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ntusoft has released a new software program for transformer 
and inductor design. The programis called Magnetics Designer. 


Magnetics Designer synthe- 
sizes and analyzes various types 
of transformers, inductors, and 
chokes. lt is optimized to design 
layer wound magnetics using cores 
that have rectangular or circular 
winding areas. Typical designs 
may include ferrite pot cores, El, 
EE, EP, cut C, and RM cores, 
toroids, and El laminated cores. 
Magnetics Designer is capable of 
designing single phase multi-wind- 
ing transformers and inductors op- 
erating from 10Hz to well above 
1MHz. Typical applications include 
60 Hz line transformers, switching 
regulator transformers, and out- 


put filter inductors. 
Continued on pg. 2 
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Figure 1, The new Magnetics Designer program makes transformer and inductor 
design easy. À demo and FAQ are available at http://www.intusoft.com. 


_ New Magnetics Design Software 


À database with thousands of cores and a wide variety of 
materials is included with the program. They are accessed 
through the Core screen shown in Figure 2. Vendors like TDK, 
Magnetics, Philips, Thomson, Micrometals, and Ferrite Inter- 
national are represented. However, the database is extensible. 


You can add your own core and material information using 
a supplied Excel spreadsheet template. 
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Figure 2, The Core Screen allows you to browse through and change the extensive 
core database included with Magnetics Designer. 


Magnetics Designer allows you to choose (Figure 1) from a 
wide variety of wire types: heavy formvar, small formvar, 
Square, double-square and Litz magnet wire (European wire 
sizes are included); pcb traces specified by the copper thick- 
ness, minimum trace width, minimum trace spacing, and PWB 
thickness and foil specified in various thicknesses. You can 
create your own wire tables with different wire sizes. For 
example, the square wire table can be edited to describe 
Virtually any rectangular wire format, or you could describe a 
table with a sparse set of gauges (ï.e., even gauges only). 


Three insulation options allow the user to specify layer insula- 
tion thickness, wrapper thickness, and end margin lengths. If 
there are multilayer foil windings, layer insulation may be 
included between each layer of foil. There is also an option to 
account for the margin required for inner layer windings to exit 
the bobbin. 
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How Magnetics Designer Works 

Tostarta transformer design, you first selectthe core family and 
material on the core screen (Figure 2). Then the operating 
frequency and electrical specifications (voltage and current per 
winding) are entered on the transformer screen (Figure 1). You 
may optionally enter one or more design constraints like tem- 
perature rise, peak flux density, and/or window fill. 


Magnetics Designer will optimize the design using a powerful 
set of state-of-the-art algorithms. The program will select the 
core geometry, turns ratio, wire sizes, and number of strands to 
achieve your design goals, all while minimizing the total power 
dissipation for a given temperature rise in the smallest possible 
core. 


Once finished, Magnetics Designer will provide a summary of 
the cores it tried to optimize, and suggestions on how to 
possibly improve the design if the design constraints were not 
met. You may browse through the core database and choose 
an alternate core geometry. lf desired, you can lock your choice 
and optimize the design around your selection. 


As shown in the transformer screen, you can manipulate over 
20 parameters such as the wire size and type, number of 
strands, gap length, insulation thickness and more. Magnetics 
Designer will predict characteristics like peak flux density, AC 
flux density, DC and high frequency AC resistance, core and 
copperlosses, leakage inductance, winding capacitance, weight, 
efficiency, temperature rise, and window fill percentage. 


Program Sports Unique Options 

Magnetics Designer has several unique features. The program 
exposes the mathematical equations and algorithms it uses to 
calculate the transformer or inductor performance (Figure 3). 
Forexample, you can access a parameter such as temperature 
rise or leakage inductance and view the equation used to 


|. variables Fr Temperature rise 
Equation 
Temperature = 710 * (Peu + Pcore) / As + Tambl 
[I NN User Defined Equation | [+] | 
Description F Datalnput 
My New Calculated = 
Temperature Data WE 
Button  : 
Units e 


Figure 3, Magnetics Designer outputs a variety of data including temperature rise, 
weight, and fillin the “User Data” Fields (left & Figure 1). The User Button dialog allows 


ou to edit the equation(s) which defines each field. 
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Figure 4, Magnetics Designer makes SPICE models for magnetics. You can 
configure and edit a schematic symbol in the IsSpice screen shown above. 


calculate its value. You can then substitute one or more 
expressions in place of the default equation. Equations can be 
set up using data from the program to derive your own perfor- 


mance merits. This feature allows Magnetics Designer to be 
customized for each application or design. 


Another unique function is the SPICE modeling capability. 
Magnetics Designer produces a Berkeley SPICE model of your 
transformer or inductor. The program also allows you to pro- 
duce and configure (Figure 4) schematic symbols. The symbols 
can be used with schematic entry packages from OrCAD 
(CAPTURES®), Protel (Schematic®) or Intusoft (SPiceNET). 


Within the general constraints described above, you can see 
that Magnetics Designer has a great deal of flexibility and can 
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Figure 5, The SPICE netlist models the core, AC/DC resistance, core loss, leakage 
inductance, and winding capacitance. Some design characteristics are also shown. 
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Figure 6, The 


general topology 


for the SPICE 


model produced 


by Magnetics 
Designer. 


handle many different winding configurations, wire types, and 
core configurations. 


How Magnetic Models Are Generated 


Figure 7, The | L___ 


wire screen 
provides 
insight into 
the expected 


losses due to |: 


harmonics 
with a plot of 
impedance 
vs. frequency 
for each 
winding. 


The complete SPICE model for a layered winding is shown in 
Figure 6. Notice that the AC resistance is approximated using 
a single L-C parallel section. The impedance is fit to the AC 
resistance curve in the neighborhood of resonance. À plot of AC 
resistance is shown in Figure 7. For all cases we have ob- 
served, the leakage inductance produces a much higher im- 
pedance than the AC resistance. This leads us to conclude that 
the effect of the AC resistance will be to damp the resonance of 
the transformer, hence the importance of making the model 
accurate in the neighborhood of the expected resonance. The 
resonant frequency using leakage inductance (L12, etc.) and 
winding capacitance represents a theoretical maximum. We 

htly by summing the previous leakage 


re 


Height  Ohms/ft 


0.2690  1.019m 
0.2400  1.286m 
02150  1.620m 
01910  2.033m 
01710  2.578m 
01580  3.250m 
0.1400  4.093m 
01260  5.156m 
01120 6.515m 


01010 8.212m 
90.10m 10.33m 
80.80m  1302m 
72.20m  1655m 
64.90m 20.70m 
58.40m  26.16m 
52.40m  3310m 


inductance and inter-winding capacitance in 
order to align our AC resistance calculation 
more closely with expected circuit operation. 


Three capacitors (Figure 8) are sufficient to 
account for the switching energy for different 
interconnection schemes. Certain turns ra- 
tios and winding directions result in 
negative capacitor values. IsSpice 
handles this properly, as long as the 
negative capacitance does not “see” 
a real valued resistor. Therefore, the 
model topologyis restricted to placing 
theinter-winding capacitances “inside” 
Figure 8, Capacitance and Ideal ofthe winding resistance and leakage 
transformer topologies. inductance models. As with most 
power electronics simulations, Gear integration works best. 
ISSPicE Users can switch to Gear integration by issuing the 
“Options parameter Method=Gear. Many other SPICE pro- 
grams, such as Pspice and Electronic Workbench, do not let 
you force Gear integration, therefore, greatly limiting their 
effectiveness for power electronics simulation. 


The ideal transformer accounts for the forward voltage transfer 
and reverse current transfer, and is modeled using a pair of 
ISSpice controlled sources (E & F). In the forward direction, the 
“primary” voltage is transferred to the secondary using a 
voltage controlled voltage source with its gain set to the turns 
ratio. In the reverse direction, a voltage controlled current 
Source assures that current is reflected from the secondary to 
the primary by the turns ratio. The voltage source after the first 
winding is used to sense current for the feedback term in the 
ideal transformer model. The DC and AC resistance equivalent 
circuits are each attached to their respective windings. 


The core is modeled as a simple resistor in parallel with the 
magnetizing inductance. The resistance is set by the core loss 
equations at the operating frequency. A linear core model is 
used. In the June 1995 Zntusoft Newsletter, We illustrated tech- 
niques for making a nonlinear saturable core model which can 
be added. The eddy current break is not sufficiently described 
in the program's core database, soits effect isn’t modeled. This 
effect can be modeled by adding a series R-C network across 
the magnetizing inductance, set to match manufacturer's data. 
This could add extra damping to the ringing observed when 
Switching devices turn off while still conducting current. Model- 
ing added effects requires considerable computational re- 
Sources, and shouldn't be done unless you depend on the 
nonlinearities for an aspect of your design. 
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Æ=  JIhe Magnetics Designer software was orig a 
ne — 1 nally developed by Analytic Artistry, a leader ANAL MIT 
intusOft in the field of nonlinear magnetics design soft- 

ware. In 1995, Intusoft acquired the rights to the “Transformer and Inductor 
Spreadsheets”, two separate DOS based programs. A large number of changes 
have been made, including combining of the DOS programs into one unified 
solution. Intusoftis nowthe sole owner and developer of the Magnetics Designer 
software. À brief list of enhancements over the DOS version is shown next: 


The core selection and transformer/inductor design algorithms have been 
changed. Magnetics Designer now automatically tries the suggestions previ- 
ously shown in the documentation. You can switch cores, families, geometries 
and materials at any time and lock your choice to play “what-if” investigations. 


IsSPice models and schematic symbols are automatically generated. 
Planar Magnetics (Printed circuit board cores and wire traces) are included. 


You can define your own core materials, geometries, and wire tables. Wire 
tables may have sparse gauges, i.e., you can have an even gauge table. 


Leakage inductance and winding capacitance values are calculated. 


The split winding feature is automatic. You can clone and move windings around 
andthe proximity fields are recalculated based on the interleaved configuration. 


AC resistance is computed for square wave harmonics (up to the 100th 
harmonic) or for a pure sine wave. 


You can set the minimum number of turns for the winding with the smallest turn 
count. You can control the accuracy of the secondary voltage outputs. 


Three thermal models are supported. The thermal model allows user entry ofthe 
convection coefficient as well as calculation of the hotspot temperature in the 
winding stack, or core, based upon heat flow through the windings or insulation 
to the core. Other models may also be implemented. 


The program supports alternate sets of units and European wire sizes. 


Magnetics Designer is compatible with Windows 3.1x, 95, or NT. 


Documentation, Pricing, and Avaïlability 

A detailed applications manual is included with the program. li provides 
information for the engineer who wants to learn more about magnetics 
design. It describes the design equations and algorithms, discusses core 
materials and geometries, and presents detailed examples. Virtually allofthe 
program’s documentation is also available on-line. 


As a Special Introductory Offer, Magnetics Designer is discounted until 
Aug. 30, 1996. (excluding North America) After that date, the price will be 
increased. Contact your local Intusoft Dealer for more information. More 
information, including an application note and a demonstration version of 
Magnetics Designer, is available on the Intusoft Web Site http://www.intusoft.com 
and CompuServe (CADD/CAM/CAE Vendor forum, Library Pi) 
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_IsSpiced Gets New BSIM3 Model 


By Farhad Aryom 


MOS IC technology has made tremendous progress in recent 
years. The desire to increase circuit speed and reduce circuit 
area has continually reduced the feature size of MOS ICs. As 
size decreases, short channel effects and high field effects, 
which did not exist in large devices, become more and more 
Significant, causing device characteristics to deviate substan- 
tially from long channel device behavior. In 1983, the Berkeley 
Short-Channel IGFET Model (BSIM) version 1 was developed. 
BSIM1 works reasonably well for devices with channel lengths 
of 1.0pm or larger. In 1989, an enhanced version, BSIM2, was 
developed. BSIM2 can fit devices down to 0.25um for both 
current, as well as small signal parameters such as output 
resistance. 


The parameter extraction strategy used for BSIM1 and BSIM2 
is relatively simple and accurate. But it introduces some inher- 
ent problems. First, a single process file does not work well 
when a wide range of device dimensions are used. Multiple 
process files are needed to obtain good accuracy over the full 
range of device dimensions. Second, the existence of redun- 
dant parameters makes statistical circuit performance 
prediction difficult, if not impossible. BSIM3v3 is the latest 
deep-submicron MOSFET model from the University of Califor- 
nia at Berkeley, and it is designed to address these problems. 


BSIMS is a physical model and is based on a coherent quasi 
two-dimensional analysis of the MOSFET device structure, 
taking into account the effects of device geometry and process 
parameters. In other words, dependencies of important geom- 
etry and process parameters such as channel length, channel 
width, gate oxide thickness, junction depth, substrate doping 
concentration, etc. are built into the model. BSIM3 allows users 
io accurately model MOSFET behavior over a wide range of 
existing technologies and predict the behavior for future tech- 
nologies. 


The BSIM3v3 (version 3) model has been extensively modified 
from previous releases and is the recommended version. Some 
of the modifications that are not found in version 2 are: 


À single I-V expression is used to model current and output 
conductance characteristics for subthreshold, strong inver- 
sion, linear, and saturation regions. This formulation 
guarantees continuity for Ids, Gds, Gm and their derivatives 
for all Vgs and Vds bias conditions. 
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There are new width dependencies for bulk charge and 
Source/drain resistance, Rds. This greatly enhances the ac- 
curacy of the model for narrow width devices. 


- _ dw and di dependencies are available for different Wdrawn 
and Ldrawn devices. This improves the models ability to fit a 
variety of W/L ratios with a single set of parameters. 


+ New capacitance equations improve the modeling of short 
and narrow geometry devices. 


+ __ New relaxation time model for characterizing the non-quasi- 
static effect of MOS circuits for improved transient behavior. 


Simulation Characteristics 

The various DC characteristics of atypical short channel BSIM3 
model are shown in Figure 9. The op-amp circuit (Figure 10) 
uses a BSIM3 model consisting of all the default parameter 
values. The frequency response for this circuit is shown below. 
À transient simulation of a one-bit comparator, along with the 
previous simulations, is also included on the Intusoft Newsletter 
subscriber’s disk, in addition to several other test circuits. 


4.00M 8.00M 


2.00M 6.00M 


@M1[GDS] 


4.00M 
) 11] 


Figure 9, The DC 

characteristics for VDS in Volts 
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TEQUENCY 


Figure 10, 
Op-amp test 
circuit using 
the new 
BSIM3 v3 
model. 


To obtain test circuits and more documentation about BSIM3V3 
and join the BSIM3 users group mailing list, you can contact the 
BSIM3 homepage at: http://rely.EECS.Berkeley.EDU:8080/ 
bsim3www/. The new BSIM 3v3 modelis included in the new 
IsSPice4 version 7.6 which is available now. 


BSIM3 Included in New 7.6 Release 


In addition to the BSIM3 model, the 7.6 release includes new 
behavioral modeling constructs and functions, a lossy trans- 
mission line with frequency dependence (skin effect/dielectric 
loss), and now R, L, C, and B elements can have expressions 
that use frequency, time and temperature. More analog and 
digital models are included (now at over 8000). Upgrades from 
previous version are available. Please contact Intusoft or your 
local dealer for information on upgrade pricing for different 
versions. 


line, MicrosimTM, the makers of the Pspice®, has stated that 
they will no longer offer DOS based products. Intusoft wel- 
comes all DOS users who still need support, and pledges not 
to force you to upgrade to Windows. We have no plans to 
discontinue our DOS based products. 


Intusoft does not force you to purchase maintenance in order to 
get full trade-up value for your existing software. Our upgrade 
policy is that upgrades are always the price difference or less! 
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In previous Intusoft Newsletters [1], we have frequently ad- 
dressed how to extend analog simulation into other domains. In 
the next few issues, we will explore how to model and simulate 
systems consisting of electrical, mechanical, and hydraulic 
elements. You might currently be using a mathematical tool like 
Mathcad or Matlab, or another non-SPICE simulator like Saber 
to analyze “mechatronic” problems, but as we’ll see, IsSpice4 
solves these problems just as efficiently with the added advan- 
tage that analog and digital electronics can be mixed in. Much 
of this work is adapted from J. Scholliers [2]. 


SPICE solves sets of nonlinear differential equations, and any 
problem that can be put into that form can be solved by SPICE. 
The added features of parametric and behavioral modeling and 


graphical input make IsSpice4 the MOST powerful and afford- 
able multitechnical tool on the market today. It should be noted 


that IsSpice4 has a full AHDL capability, but in the case of 
mechatronics, this added power is not required. 


Table 1 shows the mapping of “through” and “across” variables 
as they relate to their electrical counterparts. The energy 
conservation principle behind the relationships allows you to 
convert different systems to a single (electrical) domain. By 
translating the variables in the different domains to the electrical 
domain (voltage/current) any one or combination of the sys- 
tems can be simulated with IsSpice. Table 2 shows what 
quantities the electrical elements represent. From the user's 
point of view, once the subcircuit and symbol have been 
constructed, the electrical analog remains hidden, and the 
simulation proceeds from a graphical perspective. 


As an example of a mechatronic element, let's look at a model 
for a synchro. Synchros operate as rotating transformers, and 


CL 


2 flux rate [Wb/s =V] m.m.f. [A] 
M t S*dddt 
Mechanical force [N] velocity [m/s] F#y 
Translational F v=dx/dt 
Mechanical torque [Nm] rot. velocity [rad/s] AT 
Rotation Li © =d@/dt 
volume flow pressure [Pa] 
Fluidic F; 
[m’/s] Q 
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Table 1, 
Power 
variables 
in 
different 
energy 
domains. 


Element Mechanical Mechanical 
Fluidic 
Equation Translational Rotational 


(i.e. kg for 


Table 2, 

Electrical 

primitives voltage v=f\t) velocity velocity pressure 

in different | source (V) source (v) source ( &) source (p) 

: 1=J(L 
The units ® Si) 2) Q 
are S.I. t. dampi fluidi 
resistor (R) v=Ri - | damping (b) de PE ie 5 

() resistance 


mass, 
length) 


are used in sensitive control applications where a mechanical 
output is required. The synchro is modeled as three transform- 
ers whose gains are dependent upon position. The gain ofeach 
transformer is derived from the position of the stator with 
respect to the rotor, which is equal to K*sin(a), K*sin(o+120), 
K*sin(a-120) for the synchro (3 windings at 120 degrees). 


The rotor is modeled with a general purpose load model 
consisting of inertia, damping and Coulomb friction. This model 


* SYNCHRO - Connections: s1,s2, s3: 3 nodes in star on the stator, r1, r2: 2 nodes for the rotor, 
* w, wref: mechanical nodes. Parameters: K - gain, x0 - initial angle in radians, Jm - load inertia, b - load drag, 


* w - load friction, vO - initial mass velocity in m/s, d120 - 120 degrees in radians, K - gain 


.Subckt Synchro s1 s2 s3 ri r2 w wref {d120=2.094395)} 


xload w wref load {Jm=1f b=1f w=1f vO=0} 
xint w wref theta posrel {x0=0)} 

beqsin1 qsin1 O v=sin(v(theta)) 

begsin2 qsin2 0 v=sin(v(theta) + {d120]}) 
beqsin3 qsin3 O v=sin(v(theta) - (d120)}) 

x1 ri r2 si n qsin1 O transfov {k=1} 

x2 ri r2 s2 n qsin2 O transfov {k=1} 

x3 r1 r2 s3 n qsin3 0 transfov {k=1]} 

mnoOiT 

Ends 

* LOAD (i.e. mass with connection to ground) 
* In 101, Out 102, Gravity in m/s°2 

.Subckt Load 1 2 {Jm=1 b=in w=0 vO=0 g=0} 
xm 1 massg {Jm=1f v0=0} 

xd 12 damper {b=1f] 

xw__12fricti {w=1f} 

.Ends 

* MASSG: mass with gravity 

* Connections: 1 - X-component of velocity/force 
.Subckt Massg 1 {Jm=iu g=-9.81 vO=0] 

v 1 101 

r1101T 

r10101T 

c 101 O {Jm} ic={v0} 

ig O0 1 (J/m*g} 

Ends 

* DAMPER : Viscous damping, In 101, Ref 102 
.Subckt DAMPER 101 102 {B=1} 

r 101 102 {1/B} 

ri 10101T 

r210201T 

.Ends 

* FRICT: Constant friction force 


* 1 2 : nodes between which force has to be added 


R320101T 


.Subckt Fricti 1 2 (W=0]} 
rin 1217 

befric 1110 V= V(1,2)*{W}"0 < -.1m ? -1 : 
+V(4,2)*{W}"0 > .1m ? 1 : V(1,2)*{W}"0/.1M 

bgw 121={W}*V(111) 

.ENDS 

+ TRANSFOV: Ideal Transformator with variable ratio 
* In-pri 101 102, Outsec 201 202, Control 301 302 
.Subckt Transfov 101 102 201 202 301 302 

+ {K=1 effi=1.00 effv=1.00} 

rin 101 111 1n 

be 205 202 v={K*effv} * v(301,302) * v(101,102) 

bg 111 102 i={K} * v(301,302) * i(vm) / {effi] 

vm 205 206 

rp 101 102 1T 

rs 206 201 1n 

r3 301 302 1T 

rout 201 202 1T 

.Ends 

* POSREL: Relative Position, Input 101 (vel. of point) 
* 102(ref. vel.), Output 201 (relative position of point) 
* Parameters: x0 - initial position, K - gain 

.Subckt Posrel 101 102 201 {x0=0 K=1} 

Rin 101 102 1T 

Bw1 103 0 v={K]} * v(101,102) 

Rtime 103 104 1MEG 
Ctime 104 105 1UÙ ic={x0} 
R2 104 0 1G 

Rf 104 105 1G 

E2 0 105 1040 1G 

Rout 106 201 1N 

Eout 106 0 O0 105 1 


Figure 11, IsSPicE4 
netlist for a synchro 
model. 


.Ends 
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Figure 12, IsSPice4 simulation of 
basic synchro operation. 


can be used in other mechanical components (e.g. the rotor in 
a motor, the piston in a hydraulic cylinder). The synchro has 5 
electrical nodes (3 for stator windings, 2 for rotor windings), and 
2 mechanical nodes. The parameters for the model are the gain 
(K) of the transformer, the initial angle (x0) of the rotor, and the 
rotor parameters (inertia-Jm, damping-b, friction-w and initial 
velocity-vO). 


The rotor winding is excited by an AC reference voltage at 
frequencies ranging up to a few kHz. The magnitude of the 
voltage induced in any stator winding is proportional to the sine 
and cosine of the angle between the rotor-coil axis and the 
stator-coil axis. The voltage induced across a pair of stator 
terminals is the vector sum of the voltages across the two 
connected coils. If V sin(wt) is the reference voltage in the rotor 
of a synchro, the stator voltages are: V1 = Vsin(ot)Sin(a), V2 = 
Vsin(ot)Sin(a+120) and V3 = MAIS TO where « is 
the shaît angle. 


When an AC reference signal drives the rotor winding, the 
amplitude of the voltage induced in the stator windings varies 
according to rotor angular position. Figure 12 shows the result 
of the step change in the mechanical input along with the 
resulting synchro output. A resolver model is also included on 
the newsletter subscriber’s floppy disk. 


In future newsletters, models for a variety of building blocks will 
be developed including DC and AC motors, hydraulics, control 
elements, and sensors. This research will culminate with a 
future release of a complete mechatronic library. 


[1] "Modeling Motors”, Intusoft Newsletter #19/20, October 1989 


[2] Johan Scholliers, “Design and Simulation of Mechatronical Sys- 
tems”, VTT Automation January 1995 
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